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1 Abstract 

This paper provides results on local and global existence for a class of solu- 
tions to the Euler equations for an incompressible, inviscid fluid. By consid- 
ering a class of solutions which exhibits a characteristic growth at infinity we 
obtain an initial value problem for a nonlocal equation. We establish local 
well-posedness in all dimensions and persistence in time of these solutions 
for three and higher dimensions. We also examine a weaker class of global 
solutions. 

2 Introduction 

A fundamental question in the study of fluids concerns the possibility of 
finite time blow up of solutions to the Euler equations for a perfect, in- 
compressible fluid. It is well known that blow up cannot take place in the 
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two-dimensional case, for solutions defined over a bounded domain subject to 
Dirichlet boundary conditions (see for example the work of Wolibner [T7j and 
Ebin |7J), since smooth data in this case lead to solutions remaining smooth 
for all time. However the question remains open in higher dimensions. 

A separate class of solutions consists of those having 'stagnation-point' 
form, which attracted early attention by Weyl [TH] and Lin [13] , and provides 
a set of equations which depend on only a single spatial variable and time. 
The resulting equations, once solved, provide exact solutions to the full Euler 
equations. However, the associated growth of the full solutions in certain 
directions means that the flows possess, at best, only locally finite kinetic 
energy. Nevertheless, one may discuss such questions as finite time blow up 
for this class, and it has been shown by Stuart in [T5] that this can take place 
when the reduced equations are defined over the real line and those solutions 
decay at infinity. In this case the corresponding spatial domain for the full 
equations is M. n for n = 2 and n = 3, with the full set of solutions growing 
linearly in the other direction(s). 

The evolution of two-dimensional solutions, which can blow up on the un- 
bounded domain M 2 , therefore differs significantly from the class of globally 
defined solutions which exists for bounded subdo mains of R 2 . So the conse- 
quences of higher- dimensional stagnation-point solutions blowing up might 
be thought to result simply from their behavior at infinity rather than bear- 
ing on the question of singularity formation. This view is in some sense 
strengthened by the results of Childress et al in [3] (see also a related result 
by Cox in [5]) where solutions defined over a two-dimensional, infinite strip 
were examined. Since blow up was still found over this, smaller, semi-infinite 
domain, stagnation-point solutions defined over such domains would give the 
appearance of behaving, generally, much as those defined on the full space. 

A similar approach is implemented by Constantin in [I] to reduce the 
Euler equations, periodic in two directions, to a nonlocal Riccati equation and 
prove the blow-up in finite time by solving these equations on characteristics. 

In this paper we consider a stagnation-point class of solutions defined over 
M n which is spatially periodic in one coordinate direction. In two dimensions, 
the equations reduce to those of [3] and, although we examine slightly differ- 
ent boundary conditions, the same blow up results essentially apply. In three 
and higher dimensions we however find a 'regularizing' effect not present in 
solutions which decay to zero at infinity in the same coordinate direction, 
and this leads to the existence for all time of all such solutions, stemming 
from sufficiently smooth initial data. 
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Section [3] sets out the fundamental field equations, which in a basic sense 
date back to [13], [15]. Section @] is devoted to local well-posedness (existence, 
uniqueness and continuous dependence on initial data) of classical solutions 
to the initial value problem for the pseudo-differential equation derived in 
Section [3J We establish this result by rewriting the problem on the topolog- 
ical group T> of C 1 class diffeomorphisms as an initial value problem for an 
ordinary differential equation. Section provides apriori estimates for more 
regular classes of solutions, leading to global existence of such solutions in 
three and higher dimensions. 

In Section [6l we reconsider a class of piecewise affine solutions previously 
mentioned in [3J. These solutions are less regular than those arising in our 
existence results. It is found that they exist globally, independently of the 
underlying dimension. 

3 The n -Dimensional Equation 

Consider the n-dimensional Euler equations for an ideal, inviscid and incom- 
pressible fluid 

d t u + u.Vu + Vp = 0, 

V.u = 0, (3.1) 

where x = (sci, . . . x n ) = (xi, x'). Denoting x\ by x, u(x, x', t) represents the 
spatial velocity field of the fluid and p(x,x.',t) its pressure. We impose the 
ansatz 

u(x, x', t) = (u(x, t), — d x u(x, t)v(x', t)), (3.2) 

where the (n— l)-dimensional vector field, v, will be chosen below. As a 
consequence of (13.21) . equation (13. ip i) may be written as 

d t u + ud x u + d x p = 0, (3.3) 

together with 

d t d x u v + d x u dtv + ud 2 x u v — (d x ii) 2 v ■ V'v — V'p = 0, (3.4) 

where the primed operators refer to the variable x'. Using (13.21) and (|3.3|) . 
one sees that Vd x p = 0. Hence, differentiating (13.41) in x eliminates the 



3 



pressure term to give 

d x (d t d x u + ud 2 x u) v + d 2 x u d t v - d x ((d x u) 2 )v ■ V'v = 0. (3.5) 

Applying the V- operator to (13. 5p and using (13 .ip ii) with (13.21) . to find that 
V' ■ v = 1, shows 

d x (d t d x u + ud 2 x u) - d x ((d x u) 2 ) V'v : V'v = 0, (3.6) 

where V'v : V'v = tr(V'v) 2 = djVkdkVj (summing over j, k from 2 to n). 
For compatibility, we must choose v such that V'v : V'v is independent of 
x'. This can be done, for instance, by choosing v = — ^- x', in which case 
V'v : V'v = — and (13 .4p takes the form 

(d t d x u + ud 2 x u - ^ (d x u) 2 )v - Vp = 0. (3.7) 

In the following we examine the periodic, initial-boundary value problem, 
with boundary conditions for xGT~lR/Z, t > 0, given by 

u(0,x',t) = u(l,x',t), (3.8) 

and 

p(0,x',t)=p(l,x',t). (3.9) 

Since (13. 6 p now becomes 

d x (d t d x u + ud 2 x u) - ^ d x {{d x u) 2 )) = 0, (3.10) 

which we remark happens to be the x-derivative of a Calogero-class equation, 
(0), 

d x d t u + ud 2 u — §(d x u) = 0, 

we obtain the equation 

d t d x u + ud 2 x u - ^ (d x u) 2 = f (3.11) 
with / purely a function of time. This implies, by (13.71) . that 

V'p = ^ T x' (3.12) 

while, by (J33D, -d 2 xP = ^ (d x u) 2 + f and so Ap = --^(d x u) 2 . Finally, 
for sufficiently smooth functions u(x, t), using (13.81) and integrating (13.111) we 
have 
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f = -7TT [ (d x u) 2 dx (3.13) 
while Id, ([32]) and flU]) imply that 



dt 

Let us introduce the operator d~ l defined by 



^ I udx = 0. (3.14) 



d x 1 (f)(x,t)= (f)(y,t)dy- / <j)(y,t) dydx. 

J xo J T J xo 

We make the observation that d x and d~ l generally do not commute since 
[d x , d- l \<t> = J T (j) dx, where [P, Q\ = PQ - QP. 
Consider equation (13. lip , written in the form 

d x (d t u + ud x u) = ^ (d x u) 2 + f(t) (3.15) 

with n > 1. As a result of (13.131) and the fact that d~ l d x (f) — <f) — L <fi dx, we 
may write equation (13.151) in a nonlocal form as 

d x (d t u + ud x u) = ^ d- l d x {{d x uf) (3.16) 

and then, using the periodicity of u, we obtain 

d t u + ud x u = ^ 9- 2 ^((^ M ) 2 ). (3.17) 

Okamoto and Zhu [14] previously established local existence for (13.101) 
with u G i/ 2 , using a method introduced by Kato and Lai. Their approach 
requires showing uniqueness separately and then using uniqueness to prove 
continuous dependence on initial data. Here we instead derive a local well- 
posedness result in C 1 which follows from Picard iteration after rewriting 
the equation as an ordinary differential equation on an infinite dimensional 
Banach space. This method can also be used to prove well-posedness in 
Sobolev spaces H S (T) for s > 3/2 (see [TO] , for example, for a similar result 
for the Camassa-Holm equation). 
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4 Local Existence of Classical Solutions 



In PQ, Arnold observed that the initial value problem for the classical Euler 
equations of a perfect fluid can be stated as a geometric problem of finding 
geodesies on the group of volume preserving diffeomorphisms. Following 
this observation, Ebin and Marsden, [8], developed the functional analytic 
tools to establish sharp local well-posedness results for the Euler equations. 
This method has since been used for other equations with similar geometric 
interpretations; for example Misiolek, [12], obtained local well-posedness in 
C 1 (T) for the Camassa-Holm equation, which is the equation for geodesies 
of the H 1 metric on the Virasoro group. 

In this section we develop an appropriate analytic framework for equation 
( I3.17p . using a similar approach to prove the following theorem. 

Theorem 4.1. Suppose that n > 1. Then there exists a unique solution 

u e C ([0,T),C 1 (Tj)nC\[0,T),C°(T)), 

to the Cauchy problem for equation {3.11^ with initial data u G C 1 (T) for 
some T > 0, and the solution depends continuously on the initial data. 

Let 7 be the flow generated by u, that is, u = 7 o r y~ 1 . Then we obtain 
the equation 

7 = ^ T 9- 2 9,((9,(7o 7 - 1 )) 2 )o7 (4.1) 
from (I3.17p . Therefore it is sufficient to prove that 

^(7,7) = ^! (d~ 2 d x ((9,(7 7- 1 )) 2 )) 07 

defines a continuously differentiable vector field in a neighborhood of the 
identity on the topological group T> of C 1 class diffeomorphisms. Then The- 
orem 14.11 follows by Picard iteration over Banach spaces. 

We remark that the smooth dependence on initial data for ( 14. ip implies 
only continuous dependence on initial data for (13.171) . The map 7 — > 7 _1 
is continuous but not locally lipschitz, [8], and this prevents obtaining more 
regularity for the initial data to solution map by this method. The question 
of whether the regularity of the solution map uq — ► u(t) can be improved, 
or not, is open. It is known, for instance, that it is not possible to improve 
the regularity of this map for the Camassa-Holm equation in Sobolev spaces, 

dh- 
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In the remainder of this section, C 7 will represent generic constant depending 
only on the C l norms of 7 and 7 _1 . 

Proof of Theorem 14.11 Let us denote by P 7 the operator given by conju- 
gation 

P 7 (0):=P(0o T - 1 )o 7 

for any 7 G T> and pseudo differential operator P. Using this notation we 
write 

^(7,7) = ^t(9- 2 9,) 7 ((^) 7 7) 2 . 

Next we compute the directional derivative <9 7 P( 7i7 ) and prove that it is a 
bounded linear map. 

Note that (d~ 2 d x )f = d x 1 {f ' — f dx} is a bounded operator from C°(T) 
into C 1 (T). Furthermore we abuse the notation slightly and denote by T> the 
subgroup of orientation preserving C 1 diffeomorphisms which has index 2. 

Let s — ► 7 S be a smooth curve in D such that 7o = id and <9 s7s | s= o = W 
for W G C 1 (T). Then we have 

Tl 

d 7 F^ 7) (W) = — [ {(d £ G £ )\ £=0 o 1 + W(d x G £ ) £=0 o 1 } (4.2) 

where G £ = d x 2 d x {(d x {j o 7 e -1 )) 2 }- 

We know that d x d~ x gives the identity, hence the second summand on 
the right in (14.21) can be written as 

Wd x G £ \ £=0 o 1 =^(d x ),j) 2 ~ J^d^o^yfdx^W (4.3) 

Moreover the computation of the first summand on the right in (14.21) is re- 
duced by 

(d £ G £ )\ £=0 = d x 2 d x (d £ H £ \ £=0 ) (4.4) 

to determine d £ H £ \ £=Q where H £ = (^(7 o 7" 1 )) 2 . Here a straightforward 
computation leads to 

d £ H £ \ £=0 = - (9,(7 o 7- 1 )) 2 d x (W o 7- 1 ) - d x {(W o 7- 1 ) {d x {fi o 7 ^)) 2 } 

(4.5) 

Then, after an integration by parts, we obtain 

d £ G £ \ £=0 = - (9 X (7 o 7- 1 )) 2 (W o 7- 1 ) + (^(7 o 7- 1 )) 2 (W o ^)dx 
-d- 2 d x {(d x (^ o T l )?d x (W o 7- 1 )} 
+ (^-|){(^(7 7- 1 )) 2 (^o 7 -i)}|i 
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The last term on the left hand side of the above inequality vanishes since 7 
and W are periodic functions and 7 is an orientation preserving diffeomor- 
phism. Therefore we have 

dyF^w) = ^{ -d- 2 d x {(d x (j o 7 - 1 )) 2 a x (w o 7-1)} o 7 

+ /o 1 (^(7°7~ 1 )) 2 ^°7~ 1 ^ (4.6) 
-^/ 1 (^( 7 o 7 - 1 )) 2 ^} 

The linearity of the map W — > d 7 F^^(W) is clear. Thus we proceed 
to show that it is bounded. It is sufficient to estimate the C 1 norms of all 
three summands on the right in (14.61) . The second and third terms are both 
bounded by C 7 ||W|| c° II 711c 1 - ^or ^ ne ^ rs ^ term on the right in (14.61) . we have 

ll(^)7{((^)77) 2 (^) 7 ^}|| cl < ||(a- 2 ^) 7 {((^) 7 7) 2 (^) 7 W}|| c o 

+ ll((^)77) 2 (^)7^llco||7llc 1 (4-7) 

which is bounded by C-y 1 1 "T 1 1 1 1 Vt^ 1 1 o 1 - 

In the direction of 7, the Gateaux derivative of F is given by 

2n 

d^F^w) = -—(d-'d^m^id^w) 

and this is a bounded map since 

\\d 7 F {j>j) (W)\\ c i < C n \\(d; 2 d x ) 7 ({d x ) 7 j(d x ^W)\\ c o 

+ \\(d x )^(d x ) 1 W - [ (d^jid^W dxWcoW-yWci 

JT 

< Cn,7ll7l|ci||W||ci- 

where C nn depends only on n and C 1 norms of 7 and 7 -1 . 

In order to complete the proof of theorem 14.11 it is sufficient to show that 
F is Frechet differentiable i.e, both directional derivatives d 1 F and d 7 F are 
continuous maps. 

Continuity of (7,7) — > dyF( 7 ^(W) 

The following inequality reduces the proof of continuity of d y F^^(W) to 
estimating the two summands on the right hand side: 

\\d 7 F^ l} (W) - d 7 F (72 , 72) (W0|| cl < \\d 7 F ( ^ 7l) (W) - d 7 F (M (W)\\ c , 

+ ||(9 7 -F( 727l )(W) - d 7 F (72i72 )(W0|| c i 
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where the inequality holds up to a constant depending on n. We rewrite the 
C 1 norm that we wish to estimate to show continuity in 7 as 



1^7-^(71.7) ~ ^7-^(72 ,7) lie 1 
-2 

-1 



< 



^) 7 {(((^) 7 7i) 2 -((^) 772 ) 2 )(^) 7 iy} 
(3c (71 7 - 72 l)f W o 7 _1 (ix 



c 1 



(4.8) 



+ \W\ 



(d x (ji o 7 - 72 o j)) 2 dx 



The last two summands in (14. 8[) are bounded by C 7 ||7i — ^Hc-i || W^Hc - F° r 
the remaining term, we have 

ll(5- 2 ^) 7 {(((9,) 7 7i) 2 - ((5,) 7 7 2 ) 2 )(5,) 7 ^} Hoi 

< ll(^ 2 ^)7{(((^)77l) 2 - ((^) 7 72) 2 )(^) 7 W} Hoc 

+ 11 ((5,(7i ° 7^)) 2 " (5,(72 o 7" 1 )) 2 ) d x (W o 7- 1 ) \\co \\d xl \\ C o 

which is bounded by C 7 ||7i — 72IIC 1 ||7i + 7211c 1 II W^llc 1 - 

Our next estimate establishes continuity of 7 — > d 1 F(~ ta ){W). Note that 
it is sufficient to consider 

Wd.F^iW) - d^ id) {W)\W 

< II(5,t 2 5,) 7 {((5,) 77 ) 2 (5,) 7 w/} - d- 2 d x {{d x ^fd x w)\\ c , (4.9) 

(4.10) 

/'■'■ 

+ l|W||co 



+ 



{(5,(7 o 7 _1 )) 2 (^ o 7- 1 ) - (d x j) 2 W} dx 
{(5,(7 o 7- 1 )) 2 - (d x j) 2 } dx 



(4.11) 



where the inequality is up to a constant depending on n. After adding and 
subtracting the appropriate terms (14.101) is bounded by 

|| Vl^lloo H^-O'lloo || (^T ^" 1 )^^" 1 — ^Tllc70 + II^Tlloo ll^lloi HO'— «cZ||c70 (4-12) 

which is bounded (up to a constant C 7 ) by ||7|| 2 7 i || W^Hc 1 1|7 — ^Hc 1 - The 
term in (14. lip is estimated similarly. Hence in order to establish continuity 
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in 7 it is sufficient to bound 

II (d~ 2 d x ), {{{d^ifid^w} - d^d x ((d x jyd x w) \y 

< ||(^T 2 ^)^ {((^)^) 2 (^)^^} - 6>- 2 ^((^) 2 ^W)||c70 (4.13) 
+ 11 {(9,(7 o i^fd^W o 7- 1 )} o 79,7 - {d x ifd x W\\ c o (4.14) 

The norm in (14.141) is equal to 

HWa.7) a {(flS.7- 1 ) a o 7 -l}|| Cb (4.15) 

which is bounded by C 7 ||9 a; 7||^o||5a;W A ||c ||7 — ^llc 1 - For (14.131) it is sufficient 
to estimate 

\\(d- 2 d x ) 7 S - a- 2 ^|| c o + \\d~ 2 d x S - d- 2 d x ((d^) 2 d x W)\\ c o (4.16) 
where 

5 = 5(7,7,^0 = (d x (i °7 _1 )f 07^(^07-0 07. 
We observe that, in general, d~ 2 d x and d~ l satisfy the identities 

(d~ 2 d x ) 7 f - d~ 2 d x f (4.17) 



(CV - a- 1 / - (V Qf * f o 1 ~\x)dx^j o 7 - d- 1 Qf * 



f(x)dx 



and 



(d- i ) 1 f-d- i f=d-\f(d xl -\))- r r {x) f( V )d xl ( V , t)d V dx. (4.18) 



J X 



The second identity above is obtained changing variables in the first two 
integral terms in {d~ 1 ) 1 f — d~ 1 f . Moreover we have 



dx \J 1 ^ dx ) ° 7_ ^ X ( / f( x ) dx 

= {d x \ (J f(v)(dMv) ~ l)dv 
+(d-% ( I f(x)dx) - d- 1 ( I f{x)dx 
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if both / and d x ^f are periodic and 7 is an orientation preserving diffeomor- 
phism of T. Therefore, combining these identities, we obtain the following 
bound for the first summand in (j4.16p : 

\\S\\c4l - id\y < C 7 ||^7||^||^iy|| c o|| 7 - id\\ c i (4.19) 

The second summand in (I4.16P is equal to the C° norm of 

d~ 2 d x {d x W(d x jf {(dxj- 1 ) 3 07-I}) (4.20) 

which is bounded by C7 7 ||W||cn 11711c 1 II 7 — ^11 c 1 - Hence the continuity of 
7 — > <9F( 7 7 ) follows. 

The continuity of (7, 7) — > <9 7 F( 7i7 ) can be shown analogously. There- 
fore ^(7, 7) defines a continuously differentiable map in a neighborhood of 
(id,uo). This completes the proof of theorem 14. 11 □ 



5 Global Existence for n > 3 

In this section we investigate the persistence of solutions of the initial value 
problem for (13.171) and show that, unlike the two-dimensional case where 
solutions may blow up in finite time ([3], [5]), they persist for n > 3 in the 
appropriate function spaces. 

For the theorem below let us use the following notation, 



X n (T) 

and 

YJT) 



W 2,0O (T), n = 3 



n-l 



W 2 '^{T), n>3, 

W 1 ' 00 ^), n = 3 
W^(Ti, n>3, 



Theorem 5.1. Let n > 3 and assume that uq(x) £ X n (T). Then u(x,t) G 
C°([0,T],X„(T)) n ^([CT],^^)), for arbitrary T > 0. 

Proof Consider equation (I3.10p . expressed in the form 

d t d 2 x u + ud\u + 2=f d x ud 2 x u = 0. (5.1) 
Using the flow 7 of u (7 = u o 7) in (15. ip . we first solve for d 2 u, 

d 2 x u o j(t) = u'q exp ^— ^ — ^ J d x u o •y(s)ds^j . (5.2) 
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By the identity <9 X 7 = d x u o 7 d x j, we also have 

dxl{t) = exp (^j d x u o -y(s)ds^J . (5.3) 

Therefore 

#«o 7 (t)(0. 7 (t))£? = <. (5.4) 

Note that Theorem 14.11 implies that 7 G C 1 locally in time and it follows, 
given $3.7(0) = 1, that there exists an interval, t G [0, r(e)), over which 
< e < inf x£ j d x -y(t) < sup xe j d x -y(t) < e _1 . Equation (15. 4p then implies 
that, locally, u G X n (T), since 7 maps T diffeomorphically to itself. In 
turn, equation (13. 16[) shows that d t u G Y n (T) over the same time interval. 
(With additional assumptions on the data, further regularity can also be 
bootstrapped to higher derivatives). 

Assuming then that sufficient smoothness holds locally in time, we find 
on multiplying (15.11) by \d x u\ p ~ 2 d x u that 

d t \d 2 x u\ p + ud x \dlu\ p +p^d x u\d 2 x u\ p = 0. (5.5) 

Since equations (13.21) and ( 13. 8ft imply that both u and d x u are periodic func- 
tions of x, the same is true of d^u, by (13. lip . One therefore obtains, on 
integrating (15.51) over T, 

/ \d 2 x u\ p dx + (p2=f - 1) f d x u\d 2 x u\ p dx = 0, (5.6) 
at Jj " Jj 

71 — 1 

from which it follows that the L^^iT) norm of d 2 u is uniformly conserved in 
time for n > 3. The case n = 3 can either be considered as the limit n — > 3 
with p — > 00 in (15.61) . or directly using (15.11) which shows that d 2 u is constant 
along characteristics and hence its L°°(T) norm is uniformly conserved. 

Periodicity of u(x,t) in x implies there exists a zero for d x u, say at x — 
x (t), and so for x,x G T, 

d x u(x,t) = / d 2 u(y,t)dy. 

J XO 

For n > 3, we therefore have the estimate 

2 

|9a;u(a;,t)| < |x — xo|™ rT ||Mol|ii^i < H^oll — 

n — 3 ri — 3 
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"I 



n-1 . 

n — 3 



using Holder's inequality, and so 
If n = 3, then 

H/92 || _ II //M 

1 1 a; I 1 00 II Oil 00 

which means 

\d x u(x,t)\ <\x- XqiWu'qWoo < \ \u'q\ 
for all x G T, and so 



||5 x w||oo < |K' 



01 loo 



for all t > 0. 

Further, since u(x,t) — uq(x) has mean zero by equation (I3.14p . there 
exists x = x\{t) where u{x\,t) = Uo(xi) and, for x,x% G T, we have 



u(x, t) = u (x) + / (d y u(y, t) - u' Q (y))dy. 

It follows that 

\u(x,t)\ < ||tt ||oo + \X~ XiKII^IIoq + IKHoq) 

for all x G T, which gives the inequality 

I K| |oo < H^IU + H^oHc 1 - 
Combining the results of the previous two paragraphs shows that 

IMIc 1 < ll^ollc 2 f° r n — 3, (5.7) 

and 

\\u\lc 1 < Ikollc 1 + ||«oll2=i f° r n>3. (5.8) 

n — 3 v * 

Finally, on using the properties of the operators d~ x d x and d~ 2 d x in equa- 
tions (I3.16P and (I3.17P together with the above estimates, it is seen that 
||<9tw||oo and ||9t<? a .iz|| 0Q are majorized by a function of ||m ||c 2 for n = 3, 
while ||<9tw||oo and ||9|9 T ^||»-i are majorized by a function of II Wo 1 1 2 

for 

n-3 ' W 'n-3 

n>3. 

In both of these cases it follows that the C 1 (T) norm of u and the C°(T) 
norm of d t u remain uniformly bounded in time over any interval of local 
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existence and, by bootstrapping the arguments of Theorem 1, the solution 
can be continued, globally, in time. 

□ 

As a remark, we note here how a blow-up argument made in j3], which 
involves a nontrivial class of separable solutions to equation (13. 10|) for n = 2, 
fails to apply in the case n > 3. In particular, the possible appearance of a 
(r — t)^ 1 factor, r > 0, in the two-dimensional case no longer exists in higher 
dimensions. 

Given the solution form u(x,t) = X(x)T(t), equation (13.101) reduces to 

XX" (x) + — -X'(x)X"(x) + X(x)X"'(x) = 0, x G T, (5.9) 
n — 1 

where 

t{t) - XT(t) 2 = 0, t > 0, (5.10) 

5 — n 

and A is a constant. Multiplying equation (15.91) by \X"(x)\ r ^X"(x) now 
gives 

X\X"(x)\^ + ^(X'(x)|X"(x)|H + X(x)(\X"(x)\^y) = 0. (5.11) 
n — 1 

By using periodicity, an integration of (15.1 ip over T for n > 3 therefore shows 

/ n — 1 

A / \X"(x)\^dx = (5.12) 

and the result then follows. We note that there are in general one or more 
points of inflection in nontrivial, periodic solutions, which prevents this ar- 
gument from holding in two dimensions. 



6 Weak Solutions 

In this section, we construct a basic, piecewise differentiable class of weak 
solutions u(x,t) G C°([0, T), PC\T)) n C\[0, T), PC°(T)) to fl3~TID . which 
are found to exist for all T > 0, regardless of the underlying dimension. 

For every vector field $(x,t) G Cg°([0,T) x T x M n_1 ; W 1 ) such that 
V ■ $ = 0, and for every scalar function 0(x,t) G Cg°([0,T) x T x E" -1 ; R), 
the velocity field u(x, i) in (13. ip satisfies 
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f Q d t $ u + (V$u) • udndt + J TxRn _ 1 $(x,0) • u(x,0)dx = 0, (6.1) 

/ Rn V0-udx = O, (6.2) 

where Q = [0, T) x T x ]R n_1 . 

In terms of (I3.2p . equations (16.11) and (16. 2p reduce to 



dt<p u — <9t$' • v d x u <ix dt 
+ / 9^ii 2 - (4$' + V» ■ v^M + (V'$'v) • v^^dxdt 

JQ 

+ / 0(x, 0)u(>, 0) - $'(x, 0) • v(x, 0)d x u(x, 0) dx = 0, (6.3) 

/ <9 X u - V'0 • v x u dx = 0, (6.4) 

in which we have used the notation $ = (0, $') to distinguish the first com- 
ponent from the remaining n — 1 components of Denoting by [u] = u + — u_ 
the jump in u across any smooth surface of discontinuity, S, and considering 
test functions whose support crosses S, equation (16.2p shows that 

[u] ■ n = (6.5) 

where n = (n, n') is normal to S. Then, by (16 Ah . we have 

[u]n — [d x u]v • n' = 0, where v = r^rx'. (6.6) 

In examining weak, frontlike, piecewise continuous solutions for which [u] = 
and [d x u] ^ (see [6] ) , it follows that these discontinuities propagate so that 
n' • x' = 0. A weak formulation specific to such discontinuities may be derived 
by means of appropriate choice of test functions from (16.31) . or by observing 
that (13.1 ip may be written in conservation form as 

d t {{d x uf- n ) + d x {u{d x u) l - n ) + (n - l)(d x u)- n f = 0, x G T. (6.7) 

We will admit weak solutions, u(x,t), which satisfy the relation 

/ d t cp {d x uf- n + d x tp u{d x uf~ n -{n- l)f(p{d x u)' n dx dt (6.8) 

JQ 

+ / (p(x,0)(d x u) 1 ~ n (x,0)dx = 0, (6.9) 
Jt 
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for all <p(x,t) G C^°(£2), where £3 = [0,T) x T. Using standard Rankine- 
Hugoniot type arguments, [6], discontinuities in c^w that jump across a curve 
x — ip(t) are seen to satisfy 

(-^ + M (^,t))[(^ M ) 1 -"] = (6.10) 

and such discontinuities therefore propagate with the flow of (13. lip , i.e. ip(t) 
is a member of the characteristic family, 7 = w o 7. 

6.1 Piecewise Affine Solutions 

We begin by commenting on the general case of periodic, N— phase, piecewise 
affine solutions. Given that both d x u and d t u may be discontinuous across 
the curves x — tpi(t), 1 < i < N — 1, in order for u to remain continuous 
there we must have [u]{^{t),t) = 0, and so 4:[u](j(t) , t) = 0. As a result, 
d t [u) + ud x [u] = 0, and first derivative jumps are seen to satisfy the relations 

[d t u\ + u[d x u] = and [d x p] = 0, (6.11) 

from (13.31) . Under these conditions on u, the expression for f(t), which was 
obtained in (13.131) by integrating (13. lip for u G C 1 (T), remains unchanged: 

/ = [ {d x ufdx. (6.12) 

n-lj r 

In the special case, N = 2, which we consider here, our form of solution 
is given by the periodic extension of the function 

{xp, _ xg(o,0), 
<pp + (x - 4>) q, are(0,-0), (6.13) 

(ftp + (xp - 4>) q + (x - tjj)p, are (-0,1), 

where <p = (f>— [[</>]], tp = ip — [[ip]), with [[.]] denoting the 'integer part' of the 
argument. The functions <j)(t) and ip(t) are the representatives in [0, 1] of the 
phase curves <p(t),ip(t) e (—00,00) which start out from 0(0), -0(0) G [0,1] 
and separate those regions where d x u(x, t) periodically takes on values of pit) 
or q(t). 

Proceeding heuristically for the moment, periodicity of u(x,t) requires 
that 

M(t) = <P(t)p(t) + (^(t) - (f>(t))q(t) + (1 - m)p(t) = 0- (6.14) 
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Given the spatial periodicity in pressure (see (13. 9p ). we recall that integra- 
tion of equation (13.31) over one period showed the integral u(x, t)dx to be 
independent of time. This allows (16.131) to be used to give an expression for 
a(t). The result may be written, for instance in terms of p, and ip, as 

« + + = c (6.15) 

where we have set c = /„* u(x, 0) dx = a(0) + p(0) ~ 9(0) 0/>(O) - (f)(0)) ((f)(0) + 
-0(0) - 1). We will assume that < 0(0) < 0(0) < 1. Choosing, for 
convenience, c = 0, the 'average characteristic' must propagate with speed 
zero and we will see, consequently, that both 0(t) and 0(t) remain in [0, 1] 
for all time. The distinctions between and <f), if) and 0, will therefore not 
be made further here. 

Combining (I3.1ip . (16.121) and (16.131) now gives 

P = —,V 2 + f, q = — ~T<? 2 + / (6-16) 
n — 1 n — l 

where 

T) 

f = -Mp 2 + (0 - 0)g 2 + (i - 4>)p 2 )- (6.17) 

n — 1 

Also, by (16U01 . 

= a + 0p (6.18) 

and 

= a + 4>p+ (ip - (f))q. (6.19) 

We next verify that (16. 14f) follows from the system of equations (16.161) - 
(I6.19p . Differentiation and some simplification gives 

2 2 

M = -(^-0)(p-g)g + (V;-0)g 2 + (0+(l-^))(^- + /) + (^-0)(^- T + /) 

n — l n — l 

= -(0 - <f>)(P ~ q)q ~ ((0 + (1 - ^))P 2 + bl>- <P)q 2 ) 

= -p(q(*P - 0) + p(4> + (i - VO) 
= -pM, 

and so M{t) = Af(0)e~^o p ^ ds . In particular, taking periodic data, A/"(0) = 0, 
means that Af(t) = 0, t > 0. We may therefore use (16.141) to write (16.191) as 

(l_-0)- = _ Q;+ (l_-0)p. (6.20) 
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and again employ (16.141) to express the following phase fractions as functions 
of p and q, 

^-0 = -^-, 0+(l-V) = ^-- (6-21) 
p — q p — q 

Using these relations in (16. 1 7[) leads to 

Tt 

f = tPQ (6-22) 

n — 1 

from which (16 . 16[) reduces to the autonomous system 

p= — — -r{p 2 + npq), (6.23) 
n — 1 

q= —^—{q 2 + npq). (6.24) 
n — 1 



Subtracting (16. 18ft from (16 . 191) implies 



m - 0(t) = (^(0) - 0(0)) exp( / q(s)ds) (6.25) 

Jo 

which means that the center phase fraction does not collapse as long as 
J^q(s)ds remains bounded away from — oo. Similarly, adding (16. 18j) and 
qgHZBD gives 

<t>{t) + (1 - = (0(0) + (1 - V(0)) exp( / p(s)ds) (6.26) 

./o 

and the outer phase fraction exists as long as j^p{s)ds > — oo. Comparing 
equations (16.251) and (I6.26j) shows also that periodicity imposes the following 
requirement on p(t) and q(t) in terms of their initial phase fractions, 



(0(O) + (l-^(O))exp( / p(s)ds) + (^(O)-0(O))exp( / q(s)ds) = 1. (6.27) 

Using ( 16.161) to compute p — q next gives 

1 /"* 

p(t - g(* = (p(0) - g(0 )) exp( / p{s) + q{s) ds), (6.28) 

n-lj 

which shows p(t) — q(t) cannot change sign. By equation (16. 21 j) . p(t) and 
q(t) consequently keep their signs as long as neither phase fraction collapses. 
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This can alternatively be seen by considering a sketch of u and observing 
that p and q have opposite signs and can only vanish simultaneously. Thus, 
without loss of generality, we subsequently assume p(t) > > q(t), for at 
least some t > 0. 

In the following lemma, we solve the system (16.231) . (16.241) implicitly in 
order to show that two-phase solutions of the type (16.131) exist for all time. 

Lemma 6.1. Let (p(t),q(t)) satisfy \6. 23\) , \6-21$ with initial data p(0) > 
> q(0) (respectively, p(0) < < q(0)). Then (p(t),q(t)) exists for all 
t G (—00,00) and satisfies pit) > > q(t) (respectively, p(t) < < q(t)). 
Further, (p(t),q(t)) (0,0) as t -> ±00 and ||«(.,t)|| c i + \\d t u(.,t)\\ c o 
as t — > ±00. 

Proof Writing equations (16.231) and (I6.24p using polar variables, 
p(t) = r(t) cos9(t), q(t) = r(t) sin9(t),r(t) > 0, leads to the following, 

2 

r(t) = ^—(cos 3 (9(t))+n(cos9(t) + sin9(t))cos9(t) sin 9(t) + sin 3 (#(t))), 
n — 1 

(6.29) 

and 

9(t) = r(cos9(t) - sin 9 (t)) cos 9 (t) sin 9 (t). (6.30) 

In the case p(0) > > g(0), 9(0) G (-vr/2,0), so by (l6\30l 9(0) < and, as 
long as 9(t) G (-vr/2,0), 9 (t) < 0. 

We will show that 9(t) — > — 7r/2 as t —>■ 00 (#(£) — > as t — ► —00) by 
using ( 16.291) and ( j6.30f) . Integrating the resulting expression for gives 

r(#) = c| cos ^ sin I COS6 1 — sinflp"- 1 (6.31) 

where c > denotes a generic constant. Inserting this expression for r(9) in 
( I6.30P results in 

n+l 



9(t) = c\ COS6 1 sin^l 1 cos 9 — sin#| (cos 9 — sin 9) cos 9 sin 9 

n 2 

= — c\ cos9 sin 9\ 1 cos9 — sin 9\ n ~ 1 , (6.32) 
for 0(t) G (-tt/2,0), and so 



I cos9 - sin^l™- 1 ,„ n „. , , 

-d0 = - c t, < 0(0) < 0, c> 0. (6.33) 



'0(0) 



cos 9 sin 01 2 
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Since -^r > 1, the integral expression diverges, both as 9(t) — > —n/2 (t — > 
+00), and as — > (t — > —00). Thus and, from ( 16.311) . r(t), are 
bounded, continuous functions of time, with 6(t) G (— n/2,0). By (16.311) 
then, r(t) — > as t — > ±00 and, noting that from (16.141) . (I6.15p . 

a(t ) = - p(t) ~ g(t) (v>(t) - 0(t))(0(o) + v»(o) - 1) , 

it follows from (I6.35P and (I6.36P that 

p{t), q(t), a(t) and u(x, t) — > as t — > ±00. 
The remaing conclusions are easily obtained. □ 

Finally, we show that for c = the phases remain in the interval [0, 1] for 
all time. 

Theorem 6.1. Suppose c = and < (f)(0) < ip(0) < 1. Then the phases 
<f),ip stay in [0, 1]. In particular 

0W e [0(0), ^(0(0) + ^(0)) and V(0 e (^(0(O) + ^(O)),^(O)] 

/or all t > 0. Further, the time asymptotic behavior satisfies 

lim 4>(t) = -((p(0) + tp(0)) = lim ip(t) 

t— >oo 2 t— >oo 



Proof 

With (16.151) giving a (for c = 0), substituting into (I6.18P and (I6.20p shows 
that, by flBT26|) . 



0(t) = (l-^(t)T= |(0 + (1-V)) = i(0(O) + (l-^(O)))p(Oexp( y"p(a)cfa) 

° (6.34) 

and the individual phases therefore satisfy 

1 



0W = 0(0) + 3(0(0) + (1 - V(0)))(exp( / p(s)ds) - 1) (6.35) 
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and 



m = m - \{m + a - ^m^Mj - 1). (6.36) 

Now we examine (16.271) . Assuming q < < p, the first term is mono- 
tonically increasing in time and must converge, for < 0(0) < ^(O) < 1, to 
a positive limit. The second term is positive but monotonically decreasing, 
so may converge, as t — > oo either to a positive limit, or to zero. Formally, 
setting < J °° p(t)dt = L < oo, and — oo < q(t)dt = M < 0, equations 
fl6\27]) and flBT28|) imply 



(0(0) + (1 - m))e L + (m - 0(O))e M = 1 (6.37) 

and 

Umt^ipit) - q(t)) = (p(0) - g(0))eW. (6.38) 

If Mis finite and p(0) ^ <?(0), then (16.381) . together with the necessity for 
p{t) to approach zero as t — > oo, means that lim^oo g(t) 7^ 0. However this 
implies that M = —00, a contradiction. On the other hand, if M = —00 then, 
since L < 00, (I6.38P requires that lim^oo q(t) = 0, which is permitted. We 
conclude that, for smooth solutions pit), q(t) to exist with p(0) > > q(0), 
we require that L = J °° p(t)dt < 00 and M = J °° q(t)dt = —00. Thus, by 

eL - wrhwr (6 - 39) 

Writing (I6.35P in the form 

<P(t) = ^(0(0) - (1 - V(0)) + ^(0(0) + (1 - V(0))) exp / p( S )rfs (6.40) 







it follows that 



and similarly, 



hm0(t) = i(0(O)+^(O)) (6.41) 

hmV^) = ^(0(0) +#))). (6.42) 

By the monoticity in time of / * p(s)ds, therefore 0(t) G [0(0), |(0(O) + ?/>(0)) 
and ip{t) G (f (0(0) + ^(0)), ^(0)] for all £ > 0. 
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